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Valleytronics targets the exploitation of the additional degrees of freedom in materials where the 
energy of the carriers may assume several equal minimum values (valleys) at non-equivalent points 
of the reciprocal space. In single layers of transition metal dichalcogenides (TMDs) the lack of inver¬ 
sion symmetry, combined with a large spin-orbit interaction, leads to a conduction (valence) band 
with different spin-polarized minima (maxima) having equal energies. This offers the opportunity 
to manipulate information at the level of the charge (electrons or holes), spin (up or down) and 
crystal momentum (valley). Any implementation of these concepts, however, needs to consider the 
robustness of such degrees of freedom, which are deeply intertwined. Here we address the spin and 
valley relaxation dynamics of both electrons and holes with a combination of ultrafast optical spec¬ 
troscopy techniques, and determine the individual characteristic relaxation times of charge, spin and 
valley in a M0S2 monolayer. These results lay the foundations for understanding the mechanisms of 
spin and valley polarization loss in two-dimensional TMDs: spin/valley polarizations survive almost 
two-orders of magnitude longer for holes, where spin and valley dynamics are interlocked, than for 
electrons, where these degrees of freedom are decoupled. This may lead to novel approaches for the 
integration of materials with large spin-orbit in robust spintronic/valleytronic platforms. 


Exploiting and manipulating the spin and valley de¬ 
grees of freedom of carriers in order to process and store 
information is one of the most challenging goals of mod¬ 
ern solid-state physics, which has already resulted in the 
demonstration of several functional devices [1-6]. In this 
context, transition metal dichalcogenides (TMDs) add 
novel functionalities, due to the strong interplay between 
the spin and the crystal momentum of the carriers [7], 
and represent a promising platform to develop new spin- 
and valleytronic devices thanks to their peculiar elec¬ 
tronic structure [8] and the integrability with graphene 
technology [9, 10]. 

In a single MoS 2 layer (1 L-MoS 2 ) both the minimum 
of the conduction band (CB) and the maximum of the 
valence band (VB) are located at the K and K’ points 
of the hexagonal Brillouin zone (Fig. la) [11, 12], allowing 
for direct absorption transitions in the visible range [13]. 
The lack of inversion symmetry, combined with the 
symmetry of the Bloch wavefunctions at K, K’, lead to 
electron states with a non-vanishing projection of their 
average angular momentum (L z ) along the direction per¬ 
pendicular to the MoS 2 plane [11, 12]. In particular, at 
K, (L z ) <0 for states contributing to the CB minimum, 
while (L z ) >0 for those at the VB maximum. At K’ 
the signs are inverted, since K, K’ are related to each 
other by time-reversal conjugation [11, 12]. The valley in¬ 
dex can be thus regarded as a discrete degree of freedom 
for low-energy carriers, robust against defects, contami¬ 
nants, and low-energy phonons because of the large valley 


separation in momentum space[14], in principle enabling 
valley-based noise-resistant quantum computation [15]. 

It is important to note that (L z ) for states at the 
VB maxima is associated with the transition metal (Mo) 
d orbitals[16]. The strong spin-orbit interaction acting 
on the latter induces a AE’so =160meV energy split¬ 
ting between spin-up and spin-down states, and deter¬ 
mines opposite spin polarizations at the VB K, K’ 
(Fig. la) [11, 12]. Conversely, (L z ) at the CB minima 
is provided by a small contribution from the p atomic 
orbitals of the chalcogen (S) atoms (see Table II in 
Ref. 11), resulting in a small spin-orbit splitting (a few 
meV)[ll, 12]. States with opposite spin at the bottom of 
the CB band are thus almost degenerate, meaning that 
the spin of conduction electrons is independent of the oc¬ 
cupied valley[ll, 12]. This strong coupling between mo¬ 
mentum and spin makes single-layer TMDs particularly 
appealing when compared to other valleytronics plat¬ 
forms, such as those based on silicon[15] or graphene[17]. 
The large spin-orbit interaction, combined with optical 
excitation by circularly polarized light, allows the pho- 
togeneration of electron (e) and hole ( h ) populations in 
the CB and VB, respectively, with up to 100% spin and 
valley polarization [7, 20]. 

Thus, studying the dynamics of the relaxation of val¬ 
ley and spin polarization for both photoexcited e/h is 
crucial to understand the mechanisms leading to the loss 
of information that could be associated with these de¬ 
grees of freedom. Thanks to their peculiar electronic 
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FIG. 1. a) IL-M 0 S 2 band structure around K and K’ [11, 12]. 
The electronic states are labeled as a function of the projec¬ 
tion of their average angular momentum along the direction 
perpendicular to the plane and the spin angular momentum, 
b), Sketch of the TRCD experiment. A circularly polarized 
pump pulse, resonant with the bandgap, creates a population 
of spin-oriented e and ft in K or K’ due to valley selectivity. 
Co- and counter-circularly polarized broadband pulses probe 
the transient variation of the transmission spectrum due to 
the pump-induced phase space filling effect. The TRCD dy¬ 
namics is associated with spin-flip and valley-scattering pro¬ 
cesses (the red and blue thick arrows represent the e and ft 
spin, respectively), c), Sketch of the TRFR experiment. A 
spin polarized e and ft population is initialized by a circularly 
polarized pump pulse. The rotation angle Of of the probe 
linear polarization is related to the net value of the orbital 
momentum projection of both carriers (indicated by the red 
and blue circular arrows). Since the states at K and K’ have 
opposite orbital momentum projections, TRFR is sensitive 
exclusively to intervalley scattering processes. 


structure, TMDs offer the unprecedented opportunity of 
comparing the dynamics of VB holes, where spin and 
valley dynamics are interlocked [11, 12], with that of CB 
electrons, for which spin and valley index are completely 
decoupled[ll, 12]. This allows for the direct evaluation 
of the stabilization that the coupling between spin and 
valley provides against perturbations. 

The relaxation dynamics in TMDs was previously in¬ 
vestigated by time-resolved photoluminescence (PL) [21, 
22] or pump-probe spectroscopy [23-25]. However, due to 
the low temporal resolution in Refs. 21 and 22 and/or to 
the lack of selectivity in Refs. 23-25, these studies were 
not able to unambiguously distinguish all the scattering 
(spin and valley) relaxation channels in M 0 S 2 , and to 
decouple the e and ft dynamics. 

Here we determine the individual characteristic relax¬ 
ation times of charge, spin and valley-associated crystal 
momentum in a IL-M 0 S 2 . We find that the ft spin/valley 
relaxation time is over an order of magnitude larger than 
for e. This stems from the fact that e can exchange 
spin and valley degrees of freedom independently, while 
ft have to undergo a scattering event where spin and 
valley need to be simultaneously transferred to the lat¬ 
tice. Our results demonstrate that coupling valley and 
spin in high spin-orbit materials is a promising strategy 
towards the realization of robust spintronic/valleytronic 
platforms, where information associated with spin/valley 
degrees of freedom can thus survive for a longer time. 

To investigate the IL-M 0 S 2 spin and valley dynam¬ 
ics, we must separately access, for both e and ft, the 
spin relaxation processes occurring within the same val¬ 
ley and those involving intervalley scattering. To do so, 
we combine two complementary ultrafast optical spec¬ 
troscopy techniques: Time-Resolved Circular Dichro- 
ism (TRCD) [27] (see caption of Fig. lb for an outline 
of this technique) and Time-Resolved Faraday Rotation 
(TRFR) [19] (see caption of Fig. lb for an outline of this 
technique). We find that e relaxation, involving both 
spin-flip and inter-valley scattering, occurs on a short 
time (^200 fs), whereas ft inter-valley scattering acts 
on longer time scales^5ps. This indicates that the non¬ 
equilibrium dynamics in IL-M 0 S 2 , including the recom¬ 
bination processes, is dominated by ft scattering. 

A IL-M 0 S 2 flake (^ 10/imx30/im, see Methods) is 
mechanically exfoliated from bulk M0S2, identified and 
characterized by optical contrast, PL and Raman mea¬ 
surements, and then transferred on a fused silica sub¬ 
strate (see Methods). The use of a metal frame, as ex¬ 
plained in Methods, allows us to locate the same flake for 
both TRCD and TRFR measurements. 

In TRCD experiments the system is excited by a cir¬ 
cularly polarized pump beam, while the transient trans- 
mittivity is probed by a co- and a counter-circularly po¬ 
larized broadband probe beam. In the case of M 0 S 2 , 
the pump-induced spin and valley dynamics is given by 
the difference of the transient signals measured in the 
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two configurations: ATqb>/T( uj,t) = ATsn/T(uj,t) — 
ATqh /T(uj,t) where AT SH ( 0 h)/T is the transient varia¬ 
tion of the transmittivity measured with pump and probe 
pulses having the same (opposite) helicity (see Fig. lb 
and Methods). TRCD measurements are performed at 
77K to reduce the e-h spin scattering rate, which gener¬ 
ally increases when room temperature is approached [22]. 
The wavelength of the circularly polarized narrow-band 
pump pulse is tuned to 650nm (i.e. close to the direct 
gap optical transition measured at 77K[24]) while the 
probe pulse covers a broad wavelength range between 
500 and 700nm (see Methods). The transient optical re¬ 
sponse is dominated by the A and B excitonic transitions 
across K, centered respectively at^655 and 605nm, cor¬ 
responding to the creation of an e-h pair with e in the 
CB minimum and h either in the VB maximum (A) or 
in the split-off band (B)[28]. We observe a vanishing 
TRCD signal associated with B (Fig. 2c), while A shows 
a TRCD dynamics characterized by a bi-exponential de¬ 
cay (Fig. 2b; 655nm probe). Such behavior mirrors the 
different time-evolution of the spin-resolved bleaching of 
the direct optical transitions at K,K’, respectively, and 
can be exploited to investigate the time evolution of e/h 
spin relaxation [23, 24]. In fact, in A, bleaching can oc¬ 
cur only when the probe pulse has the same helicity as 
the pump, at least as long as photoexcited e/h main¬ 
tain their original spin and valley degrees of freedom [24]. 
The TRCD relaxation dynamics is thus determined by 
four processes (Fig. 2b): (i) h inter-valley scattering at 
the VB maximum; (ii) e inter-valley scattering and (iii) 
intra-valley spin flip at the CB minimum; (iv) e-h recom¬ 
bination. Note that spin-orbit coupling prohibits intra¬ 
valley h spin flip at the VB maximum [20]. 

Since three relaxation channels are available to e, while 
inter-valley h scattering is a more unlikely event where 
momentum and spin need to be simultaneously trans¬ 
ferred, we attribute the faster observed time scale to 
the e relaxation dynamics. The temporal dynamics of 
the TRCD signal is analyzed in terms of the rate equa¬ 
tions for the optically injected e/h spin densities, re¬ 
spectively e^)(£) and h^^(t). In particular A T/T(t) is 
proportional to the product between unoccupied states 
in the CB [ucb — e i(t)\ an( 4 occupied states in the VB 
[nvB — h^(t)] where ticb/^vb are two parameters ac¬ 
counting for the available states in the CB/VB. The 
fitting function used to simulate the TRCD trace in 
Fig.2b is: AT/Tcr>(t) oc {[ncB — ej,(t)] [nyB — Aq(t)]} — 
{[^cb — e^(£)] I n VB — hq(£)]} (see Methods). It is impor¬ 
tant to stress that the above-mentioned mechanisms af¬ 
fecting the e/h dynamics are indistinguishable from the 
point of view of TRCD and can be described by an ef¬ 
fective relaxation time r e = (l/r e _i n t ra + 1 /^e-inter + 
1/t 6 — rec) ^ lOOfs for e and Tfo — (I/t/i— inter T 
1/Ri— rec ) —1 ~ 4ps for A, T e _i n t r a being the time con¬ 
stant of the e intra-valley (spin flip) relaxation, r e _ i nte r 
{jh- inter) the characteristic time of the inter-valley e ( h ) 
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FIG. 2. a) AT/T spectra at 77K for co- and counter- circu¬ 
larly polarized probe pulses (red and blue traces, respectively) 
at different delay times. The circularly polarized narrow-band 
pump pulse is tuned to 650nm, while the probe pulse covers a 
broad energy range between 500 and 700nm. b), Black dots: 
difference between the AT/T signal at 655nm (peak of the 
A exciton) measured by co- and counter-circularly polarized 
probe pulses and normalized to unity. Red solid curve: fit 
with the rate equations describing the temporal dynamics of 
e tU-)W an d hf(\.)(t). Inset: sketch of the processes causing 
the time-evolution of the TRCD signal, c) Comparison be¬ 
tween the difference between the spectral traces at r = 0.1 ps 
(black dots) and the spectrum of the probe pulse used in 
TRFR 


scattering, and r e (h) -rec the e (h) radiative and non- 
radiative recombination time. The value of is in agree¬ 
ment with the spin/valley decay time (^ 4 ps) measured 
by time-resolved PL [22] , which sets a lower limit for the 
e/h recombination time. This reveals that recombina¬ 
tion does not affect the e dynamics, which is dominated 
by fast intra- and inter-valley scattering, and that polar- 
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ized PL is given by the radiative recombination between 
unpolarized e and almost fully-polarized h. 

We further investigate the inter- and intra-valley dy¬ 
namics governing the non-equilibrium response of 1L- 
M 0 S 2 with TRFR, where a left (right)-circularly polar¬ 
ized pump pulse, resonant with the A exciton, creates 
a spin- and valley-polarized density of e/h with non¬ 
vanishing and well-defined orbital angular momentum. 
A delayed linearly polarized probe beam is transmitted 
through the sample, subsequently separated through a 
Wollaston prism into two orthogonal linearly-polarized 
beams and equally distributed on two balanced photo¬ 
diodes. The transient dynamics associated with the or¬ 
bital momentum of the photoexcited charges is probed 
by measuring the rotation angle Op of a linearly polar¬ 
ized probe pulse, which manifests itself as a difference 
between the signals measured by the two photodiodes 
(Fig.lc). Because of the direct gap transition at the two 
degenerate valleys, the TRFR signal in M 0 S 2 is affected 
both by Pauli blocking, inhibiting the absorption of light 
having the same circular polarization of the pump, and 
by helicity-dependent light scattering from the photoex¬ 
cited e/h populations. The first phenomenon is similar to 
that governing TRCD: e/h populations with well-defined 
spin and valley degrees of freedom reduce the absorp¬ 
tion of probe light with the same circular polarization as 
the pump. Therefore, the imaginary part of the tran¬ 
sient refractive index assumes different values for oppo¬ 
site light helicities. Since the real and imaginary part of 
the refractive index are related by the Kramers-Kronig 
transformations [29], the system displays asymmetric val¬ 
ues of the phase accumulated by left- and right-circularly 
polarized light crossing the sample, resulting in the ro¬ 
tation of the polarization plane of the linearly polarized 
probe light. This process, however, is active only when 
the probe light is resonant with electronic transitions 
across the gap. In order to rule out this effect, we per¬ 
form a two-color TRFR experiment in which the probe 
energy is deliberately tuned below the absorption edge 
(A P robe = 700nm) as shown in Fig. 2c. 

Fig. 3a plots the dynamics of the Faraday rotation, 
Op(t), measured at 77 K by exciting spin polarized 
charges in either the K or the K’ valley with a circu¬ 
larly polarized pump pulse. The change of sign of the 
signal with excitation helicity is due to the valley selec¬ 
tivity of the near gap resonant transitions. As expected, 
TRFR performed with a linearly polarized pump pulse 
(gray curve) gives rise to a negligible signal because both 
valleys are equally populated. 

As far as electric dipole transitions are concerned, pho¬ 
tons do not couple with the spin of the charges but only 
with their orbital degree of freedom [19]. Thus, the Fara¬ 
day rotation associated with the scattering of light by the 
photoexcited charges can only derive from an unbalanced 
distribution of their orbital angular momentum projec¬ 
tions. As shown in Fig. 3b, the orbital momentum of both 




FIG. 3. a) TRFR traces measured by exciting the sample 
with left and right circularly-polarized pump pulses (red and 
blue curves, respectively) and a linearly polarized pump pulse 
(gray curve). Excitation is centered at 650nm and the probe 
at 700nm. All the traces are measured at 77K. The inset 
shows the temporal evolution of 0p(t) in semi logarithmic 
scale in order to emphasize the bi-exponential character of 
the decay rate. The black solid line represents the analytical 
fit of the data to the rate equations. The decay constants 
Te-inter and Th -inter are related to the intervalley scattering 
of e and h, respectively. These processes are accompanied by 
an exchange of crystal momentum and are sketched in b). 

e and h is associated with the valley degree of freedom, 
therefore the TRFR experiment is sensitive exclusively 
to inter-valley and recombination dynamics, and it can¬ 
not be influenced by intra-valley transitions in the CB 
minimum. 

We now analyze the TRFR dynamics using the rate 
equations for the photoexcited e K ( K '){t) and h K ^ K ')(t) 
densities in the K(K’) valley (see Methods). The TRFR 
trace is then given by Spit) ~ [A e ex(t) + AhhK(t)\ — 
[.A e eK'{t ) + AhhK'(t)\ where A e ^) are two fitting coeffi¬ 
cients, taking into account the different weights of spin- 
oriented e and h (see Methods). As shown in the inset of 
Fig. 3a, Op(t) relaxes with two different time scales. We 
attribute the faster one to e inter-valley scattering, with 
a decay time r e _i nte r ^ 200 fs, while we assign the slower 
one to h relaxation with ~ 4 ps, in excellent agreement 
with the TRCD results. With a reliable estimation for 
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t e—inter and considering that the recombination process 
are markedly slower than all the other e relaxation pro¬ 
cesses (i.e. 1 /t e—rec ^ 0), we can now extrapolate from 
the TRCD data the characteristic time constant for intra¬ 
valley e spin-flip transitions, obtaining r e _i ntr a 00 200 fs 
at 77 K, comparable to r e _ i nte r- 

We then study the temperature dependence of the 
TRFR dynamics (Figs.4a,b,c). We ascribe the observed 
decrease of r e _i nte r and with temperature to the in¬ 
crease of the momentum scattering rate by phonons. The 
pump fluence dependence of 0p(t) is also investigated in 
Figs.4d,e,f. Both decay constants (r e _i nte r and 77 J de¬ 
crease exponentially for an increasing density of photoin- 
duced charges, contrary to the case of lL-WSe 2 , where 
it was shown that exciton-exciton interaction does not 
affect the spin relaxation time [25]. We attribute this 
behavior to the exchange interaction acting on the e/ft 
populations, which increases with pump fluence. We ex¬ 
trapolate the value of r ^-inter from for E pump —> 0, 
obtaining r h - int er,o ~ 5 ± 1 ps. 

In conclusion, we determined the inter-valley and intra¬ 
valley scattering dynamics of optically injected spin- 
oriented charges in IL-M 0 S 2 by combining two comple¬ 
mentary techniques, TRCD and TRFR. We found that 
the characteristic time scales of spin-flip and inter-valley 
scattering for electrons are both~ 200fs, whereas hole 
inter-valley scattering acts on a time scales 5ps. The 
large stabilization of spin/valley degrees of freedom re¬ 
sulting from the coupling between spin and valley dy¬ 
namics paves the way to the integration of materials with 
large spin-orbit in robust spintronic/valleytronic plat¬ 
forms. Furthermore, the procedure described here is of 
general validity and can be applied to study e/ft intra 
and intervalley scattering processes in any other mono- 
layer TMD. 
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Ricerca grant No. RBFR12SW0J of the Italian Ministry 
of Education, University and Research and SEARCH-IV 
No. 2013-0623, grant of the Fondazione Cariplo. 

METHODS 
Sample Preparation 

M 0 S 2 hakes are produced by micromechanical cleavage 
of bulk M 0 S 2 crystals (Structure Probe Inc.-SPI, natural 
molybdenite) onto Si wafers covered with 285nm Si02- 
IL-M 0 S 2 hakes (Fig. Mia) are then identified by optical 
contrast, PL and Raman spectroscopy [30, 31]. Fig.Mlb 
shows a representative Raman spectrum of the flake used 
for our experiment, acquired at 514.5nm with a Renishaw 


microspectrometer. The position of the two main Ra¬ 
man peaks is^ 385 and^ 404cm -1 , indicating that this 
is IL-M 0 S 2 [30-32] . Fig.Mlc plots the PL spectrum of our 
flake measured at room temperature for 514.5nm excita¬ 
tion. This consists of two bands at^ 1.85 and~ 1.98eV, 
consistent with the A and B excitons in 1 L-MoS 2[33]. 
The selected IL-M 0 S 2 hake is then moved by a wet- 
transfer technique based on a sacrificial layer of poly¬ 
methyl-methacrylate (PMMA)[34, 35]. The polymer is 
deposited onto the flake by spin coating, followed by im¬ 
mersion in de-ionized water. Water intercalation at the 
PMMA-Si02 interface detaches the polymer film[34, 35], 
with the IL-M 0 S 2 hake attached on it. This is then 
moved onto a 100 fi m thick fused silica substrate and 
left to dry for a few hours. The fused silica substrate 
is particularly suitable for TRFR experiments since it 
is an isotropic material, without crystal orientation, and 
it does not exhibit birefringence [36]. PMMA is then re¬ 
moved by acetone, and the flake is released onto the fused 
silica [34, 35] . A metal frame is then fabricated around the 
selected IL-M 0 S 2 by photolithography, followed by ther¬ 
mal evaporation of 2nm Cr and lOOnm Au films. This 
ensures the same flake can be easily found for TRCD and 
TRFR measurements. To ascertain that no damage or 
changes are induced on the selected flake by the trans¬ 
fer process, this is further characterized after transfer on 
fused silica and after the metal frame is defined, but be¬ 
fore TRCD and TRFR. Fig. Ml shows the selected 1L- 
M 0 S 2 transferred on glass and its Raman and PL spec¬ 
tra. These show no significant changes with respect to 
the initial sample prior to transfer. 

Time-resolved Faraday rotation 

A regeneratively amplified mode locked Ti:Sapphire 
laser, providing 150 fs, 500/iJ pulses at 780nm and 1kHz 
repetition rate drives two optical parametric amplifiers 
(OPA). The output of the first OPA, used to inject spin 
polarized carriers in the sample, is tuned to 650nm, i.e. 
quasi resonant to the bandgap wavelength, and it is cir¬ 
cularly polarized by a quarter wave plate. The probe 
pulse, provided by the second OPA, is linearly polarized 
and it is centered at 700nm. Both pump and probe pulses 
have a spectral bandwidth of 10 nm. The temporal de¬ 
lay between the two pulses is adjusted by a mechanical 
delay stage. At the sample position the pump (probe) di¬ 
ameter is 50(20)/im. We estimate a temporal resolution 
Ttrfr of the TRFR setups 70fs. The sample is posi¬ 
tioned in a liquid nitrogen cryostat and the temperature 
is checked by a thermocouple placed in proximity to the 
sample. The transmitted probe pulse passes through a 
Wollaston prism and it is focused on a couple of balanced 
photodiodes. The Wollaston prism is rotated in order to 
equalize the probe intensities on the two photodiodes. 
The pump induced imbalance of the signal measured by 
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FIG. 4. a) Normalized TRFR dynamics measured at different temperatures in the 77-300K range. b,c) Temperature 
dependence of the decay constants r e -inter and th- inter- d), Normalized TRFR dynamics measured at different pump fluences 
in the 1.4-20nJ range. e,f) Fluence dependence of r e - inter and rn- inter- th- inter ~ 5 d= 1 ps is extrapolated from Fig.4f for 
-Fpump —>-0, by deconvoluting with respect to the carrier lifetime 


the two photodiodes is registered by a lock-in amplifier 
which is locked to the modulation frequency of the pump 
beam (i.e. 500Hz). 


Time-resolved circular dichroism 


The pump beam is generated by an OPA at 650nm and 
it is circularly polarized. The probe is a white light pulse, 
ranging between 500 and 730nm, generated by focusing 
the fundamental pulse in a 2mm-thick sapphire plate. 
The probe pulse is circularly polarized by a broadband 
quarter wave plate. We estimate the temporal resolution 
of the TRCD setup to be^ 200fs for the entire spec¬ 
tral window of the probe beam. A fast analog-to-digital 
conversion card with 16-bit resolution enables single-shot 
recording of the probe spectrum at the full 1kHz repe¬ 
tition rate. By recording the transmitted spectra with 
(T on ) and without (T 0 ff) the pump excitation at differ¬ 
ent delay times one obtains a two-dimensional differen¬ 
tial transmission (AT/T) map as as a function of probe 
frequency uj and delay r: 


A T(u,r) 


T oxi {u3,t) - Tpff (u) 


Rate equations for spin oriented electron and hole 
densities 


TRCD results can be interpreted in terms of the time 
evolution of the spin oriented e and h densities, which 
are optically injected in IL-M 0 S 2 through optical orien¬ 
tation. Let e^)(£) be the e densities generated at K, 
K’ and h^^(t) the h densities at K, K’. We describe 
the temporal evolution of the spin polarized e/h densi¬ 
ties, following an excitation with left-circularly polarized 
light, with a set of four coupled rate equations, similar 
to those used to model the M 0 S 2 exciton relaxation rate 
in Refs. 20 and 37: 


dei(t) 

dt 


= S(t)~ (ei(t) - e t (f)) 
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de^(t) 
dt 
dh^(t) 


dt 
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The gaussian function S (t) represents the laser source 
term which accounts for the temporal resolution of the 
TRCD setup. The transient depletion of the ej. and //-j 
carriers (Eqs.Ml, M3), injected by the pump in a val- 
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FIG. Ml. a) Selected IL-M 0 S 2 flake on SiCb- b) Raman 
spectrum of the flake used for our experiments for 514.5nm 
excitation, c) PL spectrum of the selected flake measured for 
514.5 nm excitation measured at room temperature. 


FIG. M2, a) Selected IL-M 0 S 2 flake after transfer on fused 
silica, b) Raman spectrum of the flake in (a) measured for 
514.5 nm excitation, c) PL spectrum of the flake in (a) for 
514.5 nm excitation measured at room temperature. 


ley, and the consequent filling of the reversed spin po¬ 
larized carriers and (Eqs.M2, M4) are regulated 
by spin-flip processes. For e these processes can occur 
either in the same valley, or by scattering in the oppo¬ 
site valley. The former scattering process is regulated 
by the time constant r e _i ntra , while the latter depends 
also on r e _i n ter- For h spin and valley relaxation coin¬ 
cide, and spin reversal processes are accompanied by an 
exchange of linear momentum (r^- i nte r)- V^Ti-rec is the 
h recombination rate, while the e recombination rate is 
one order of magnitude slower than the intra- and inter¬ 
valley scattering rates, thus it can be safely neglected in 
Eqs.Ml, M2 (l/r e _ rec ~ 0). Once the Eqs.Ml-M4 are 
solved, the TRCD curves are generated by the product 
of the unoccupied states at the top of the valence band 
and unoccupied states at the bottom of the conduction 


band: 


—^~ {[^cb - ei(t)\ [n VB - h^(t)]} + 

- 1 J CD 

- {[^cb - [ny b - hi(t)]} , (M5) 

where tiqb and ny B are free fitting parameters account¬ 
ing for the available states in CB and VB. 

Considering TRFR, the set of four rate equations has 
to be slightly modified in order to take into account the 
sensitivity of the measurement with respect to the pho¬ 
togenerated densities of e/h in K and K’. Let e K ( K ')(t) 
be the e densities generated at K, K’ and hjc(K')(t) the 
h densities at K, K’. In this case, considering an exci¬ 
tation with left-circularly polarized light, the set of four 

























coupled rate equations can be written as: 


de K {t) _ $ _ e K (t) - e K '(t) 

dt Us—inter 

dew jt) _ e K {t) - e K '{t) 

dt 7e —inter 

dh K (t ) _ g ^ _ h K {t ) - h K '{t) h K {t 
dt r h 

—inter Ui—re 

dh K '(t) _ h K (t) - h K '(t) h K '(t ) 

M ^ ^ 


Ui— inter 


(M6) 

(M7) 

(M8) 

(M9) 


where the full-width half-maximum Ftrfr of the gaus- 
sian function takes into account the temporal resolution 
of the TRFR setup. The TRFR is only sensitive to the 
e/h intervalley scattering time, which affects the unbal¬ 
ance of orbital angular momentum, optically injected by 
circularly polarized light so that the TRFR fit function 
can be written as: 


# f (£) ^ [A e eK(t) + AhhK(t)\ — [A e eK'{t ) + AhhK f (t)] , 

(M10) 

where A e (/q are two free fitting parameters whose ratio 
is varied in order to obtain the best fit of the data. The 
simultaneous fit of the TRCD and TRFR traces with 
Eqs.M5, M10 enables us to determine the time constants 
for intra- and inter-valley scattering for both electrons 
and holes. 
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